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Using an exact approach, we study the dynamics of entanglement between two qubits 
coupled to independent reservoirs and between the two, initially disentangled, reservoirs. 
We also describe the transfer of bipartite entanglement from the two-qubits to their 
respective reservoirs focussing on the case of two atoms inside two different leaky cavities 
with a specific attention to the role of the detuning. We present a scheme to prepare 
the cavity fields in a maximally entangled state, without direct interaction between 
the cavities, by exploiting the initial qubits entanglement. We discuss a deterministic 
protocol, working in the presence of cavity losses, for the generation of a W-state of one 
qubit and two cavity fields and we describe a probabilistic scheme to entangle one of the 
atoms with the reservoir (cavity field) of the other atom. 
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1. Introduction 

The decoherencc of entangled systems has been the subject of an intense research 
activity in the last few years, triggered by the theoretical prediction of the phe- 
nomenon of entanglement sudden death (ESD) LU, experimentally demonstrated 
with entangled photon pairs ^ and atomic ensembles ™ Subsequent theoretical pa- 
per have extended the original work in various directions, by studying the dynamics 
of entangled systems coupled to either the same ^ or different reservoirs, includ- 
ing non-markovian effects ^ experimentally tested in 121 by examiningthe effect of 
finite temperature reservoirs ^ by including count er-rot ating terms ^ and by ex- 
ploring the role of diversity between the subsystems I ' \ Furthermore, it has been 
understood that the sudden death phenomenon is, in fact, rather a transfer, as 
the environ men tal degrees of freedom become entangled while the two subsystems 
disentangle L^ (entanglement sudden birth, ESB). 
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The study of reservoirs entanglement in relationship to the phenomenon of en- 
tanglement sudden de ath has been also performed for two-level-systems interacting 
with two ideal cavities ' ' \ Further investigations have been devoted to the exis- 
tence of entanglement invariant LLal 

The first aim of this paper is to study the sudden birth phenomenon in the case 
of non-markovian reservoirs, in order to elucidate the role of entanglement memory 
on entanglement exchange. In doing this, we will show that the two cases of «)a 
spectrally flat, marko yian en viron ment, and of ii)& monochromatic environment, 
first studied in Refs. ^^ and LL^, respectively, are just two manifestation of the 
sam e ph enomenon, as both of them can be described in a unique fashion (see also 
Ref. LL21). Indeed, the Hamiltonian model that we discuss continuously interpolates 
between these two extrema depending on the spectral width of the environmental 
projected density of states. 

More explicitly, we present an exact approach to the dynamics of two qubits 
interacting with independent structured reservoirs having Lorentzian spectral dis- 
tributions. This system describes, e.g., two atoms inside two leaky cavities, whose 
quality factors are related to the widths of the Lorentzians, and provides a quite 
reali stic description of the case of two ions trapped into two cavities, as discussed in 
Ref. L^. The solution of this model reduces to the Markovian one in the bad cavity 
limit and to the Jaynes-Cummings one in the ideal cavity limit. 

The exact solution of the total closed system allows to describe the zero- 
temperature reservoirs as effective two-state systems, therefore the entanglement 
between the qubits, between each qubit and its own reservoir and between the 
two reservoirs can be characterized by using the concurrence 1^4 In this frame- 
work we show how, for the general non-Markovian case, the disappearance of the 
initial atomic entanglement in a finite time (ESD) is accompanied by the sudden 
appearance of entanglement in the reservoirs (ESB), as already demonstrated in 
the Markovian case, and discuss with particular attention the role played by the 
detuning of the atoms with respect to the cavity fields. Moreover we prove that the 
qubit-reservoir tangle is an entanglement invariant. 

Finally, we describe a deterministic scheme to manipulate the system in the 
presence of losses, in order to generate a W-state of one atom and two cavity 
fields. Starting from this state, and for certain values of the system parameters, 
the subsequent measurement of the logical state of the qubit projects the cavity 
fields in a maximally entangled state. Moreover we demonstrate how one of the two 
atoms can be maximally entangled with the field inside the other cavity. Similarly 
to the case of entanglement swapping, we show that quantum systems that have 
never interacted directly, namely one atom and the field inside a distant cavity, 
can be entangled if some entanglement is present in the initial two-atom state. As 
already mentioned, these procedures work independently of the amount of non- 
markovianity of the environment; that is, independently of the quality factors of 
the two cavities. 

The paper is structured as follows. In Seed] the analytic solution of the model 



July 4, 2011 0:13 WSPC/INSTRUCTION FILE SREntangv2 



Non-Markovian dynamics of system-reservoir entanglement 3 

is presented together with and its Markovian and ideal cavity limits; Sec. |3] shows 
the main features of the entanglement dynamics; in Sec.|4]the scheme for entangled 
state generation is reported, while a summary and some final remarks are given in 
Seed 

2. Exact dynamics 

Consider a pair of two-level systems initially entangled to some extent and inter- 
acting with two different environments; a situation realized, e.g., for two initially 
excited (and entangled) atoms spontaneously emitting into two different cavities. 
Each of the two atom-cavity systems is described by Hamiltonian of the kind 

H^=a+J29k a h a e i5k " t + h.c. (1) 

k a 

This Hamiltonian is written in the interaction picture with respect to the free atom 
and field terms. Here a = 1, 2 distinguishes between the two atom-cavity system, 
<Ta are the rising and lowering operator for the atom a spontaneously emitting a 
photon into the k a modes, for which bk a are the annihilation operators. Finally the 
detuning is defined as the difference between the atomic and the mode frequencies 

Sk a =W a - LOk a ■ 

Each of the two sets of cavity modes is specified by its spectral density and time 
correlation function: 

U{t) = Y, \9kJ 2 e iS *° T = I ' dojJ a (oj)e i5 ^ r . 

k a 

We will assume that each cavity has a finite photon lifetime and describe this by 
assuming Lorentzian broadened spectral densities 

O 2 \ 

j (, A — _ a ^2 

-K (LO - UJ C J 2 + \ Z a 

where cj Cq is the central frequency of cavity a, while A^ 1 is the cavity life-time. 
For future reference, we also define the detuning of the atomic transition from the 
frequency of the main cavity mode: S a = ui a — uj Ca . As one can evaluate using the 
residue method, such density of states leads to a correlation function of the form 

-00 * 

from which we deduce that j- is the correlation time for the a-th reservoir. 

For different values of the parameter A Q , this model interpolates between an 
ideal Jaynes-Cummings model with vacuum Rabi frequency fl a (which is obtained 
for X a — > 0) and a Markovian environment giving rise to an exponential decay of 
the atomic excited level with a rate j a = 2Q,^/X a , which is obtained for a very 
large A Q . For intermediate values of the lifetime the model is known to describe 
finite memory effects (with a memory time given by the photon lifetime within the 
cavity) . 
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Assuming that the two cavities are initially empty and that the two atoms are 
prepared in an entangled state, we have 



1^,0) 




\9192) 



,%<p 



1 



eie 2 ) 



|0 fc J 



(2) 



a=l,2 



where s G [—1, 1] is a separability parameter, related to the initial tangle as 
Ti,2(0) = 1 — s 2 . Notice the the sign of s is irrelevant for the entanglement, as it only 
determines the initial population inversion of the two atoms (that is, (cr^(0)) = s). 
The time evolution generated by H^ l > gives rise to 



hM) 



1-5 



3i22> (g) |0O + e 



i£V>®l02> 



where 



b a ) = £ a (t) \e a ) <E) \0 ka ) + F a (t) \g a ) (g> |photonJ , 



(3) 



(4) 



with |photon Q ) being the (normalized) state describing the photon emitted in the 
modes of the a-th reservoir, whose explicit expression is not relevant. 
For reservoirs with the spectral density given above, one has: 



£ a {t) = e 



(\ a -iS a )t 



,,w a t X a - iS a . ,,w a t 

cosh — ) + smh- -) 

2 w„ 2 



with 



\F a {t)\ = y/l-\£ a (t)\*, 

^\l - 2is a x a - mi 



(5) 
(6) 



«'o 



where R a = ^/n 2 a + 51/4. 

2.1. Limiting cases 

On resonance, S a = 0, it is easy to evaluate the limiting expressions for the ampli- 
tude £ . In the Markovian limit, it takes the approxi mat e form £ a (t) = exp{ — j a t/2} 
and the problem reduces to the one described in L^M On the other hand, in the 
Jaynes-Cummings limit one simply has £ a (t) = cosil Q i. 

Some more details about these two limiting cases are reported in the following 
subsections. 



2.1.1. Markovian limit 

Taking the limit A Q — > 00, one has w a ~ A Q (1 — p 2 -), and furthermore, the hyper- 
bolic functions can safely be approximated by the positive exponentials so that 



exp{-^} = e-^ 1 ' 2 



(7) 
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This could have been seen directly in the spectral density: 

-,. t / \ t "fa "a ^a "fa 

Inn J a (u) = hm — — = — , 

A„-s-oo A a -s-oo 27T (W — 0J Ca ) ~r ^a ^ 7r 

which implies that / a (r) = 7 Q <5(r) . An approximately flat spectral density leads 
to a short memory (markovian) reservoir. The approximate flatness of J a (uj) also 
leads to the fact that the detuning is irrelevant as long as it is smaller than the 
reservoir band-width. Indeed, by including the leading term in 6 a one has 

£ a ~ cxp{- f A + ^) t} (8) 



4A C 



2.1.2. Jaynes-Cummings limit 

The ideal cavity limit is obtained for X a — > 0. By just dropping all of the terms 
containing A Q , one has 

i5 a 



f 



iS a t/2 



cos R a t 



2R C 



sin R a t 



(9) 



Including the leading correction in A, the expression for £ becomes quite cubersomc; 
while it considerably simplifies in the resonant case (S a = 0) to become 

XI 



-X a t/2 



cos Q a 



Sila 



t. 



3. Entanglement dynamics 

The general expression for the concurrence as a function of time is 

Ci, 2 (t) = max JO, |£i£ 2 | Wl - « 2 - (1 + *) l^i-^l] } , 
and the entanglement sudden death is found to occur at a time such that 



(10) 



l-FiW-Fatol 



1 - s 



1 



Since |.F| < 1, this implies that s > is required, which simply indicates that, for 
ESD to occur, the component of the initial state taking part to the decay should 
be larger than the other, stable one involving the ground states. As noted in l^l, 
however, due to the memory of the reservoirs, the death is (partly) reversible and 
the entanglement can revive (and then die again) after some time. In the limit of 
perfect cavities, A — > 0, the entanglement dynamics is truly periodic with a sequence 
of dark periods followed by perfect entanglement revivals. 

Since the state given above is in the Schmidt form, the reservoir behave as 
effective two level systems. This implies that, at any given time, one can evaluate 
the entanglement between the reservoirs by employing the concurrence again. It is 
easy to show that 



^i"i,r2 \t) 



{o,| 



^1^2 
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(ii) 
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It is clear from the form of the initial state that there is no entanglement between 
n and T2 at the beginning of the time evolution, but such an entanglement builds 
up after some time, and a "sudden birth" is found to occur. 

It also very easy to evaluate the entanglement between qubit a and its reservoir. 
Since the reservoirs are effective two level systems, the concurrence can be used once 
more, and one hnds: 



C' 



a 



a)\£ a F a \. (12) 

Furthermore, one can use the tangle to evaluate the entanglement between the two 
(qubit+cavity) systems. It is easy to verify that 



-s)\£ a \ 2 /2 {l + s)£ a T*j2 
s)£* a T a /2 {I + s)\T a \ 2 /2 
(1 



Pa.r 



o (1 
(1- 

Vo 



\ 



0/2/ 



so that the tangle 



r (l®ri)/(2®r 2 ) 



2(l-Tr{p? iri })=l 



(13) 



which is equal to the square of the initial concurrence between the qubits. Thus, 
we have found a simple entanglement invariant: the full Hamiltonian is bi- local 
with respect to this bi-partition of the overall system, and therefore this type of 
entanglement doesn't change in time. Notice that this is true irresp ective of the 
form off, so that it holds true for the Jaynes Cummings model! ' 'as well as for 
the markovian decay L^M. 

While for entangled states with a single excita tion som e conservation law has 
been found for the squares of the concurrences, see I I If, for the two-excitation- 
case we consider here, only the inequality 



C 



<T< 



< C 12 + C rir2 + C 1>ra -I- 2ri ^ 7"(lig>n)/(2<8>r3) 

has been verified, while a m ore general relation (not directly involving the concur- 
rences) has been reported in^^. As a first example of the time evolution described 




/.\ . ./v. ....... -'. nt 




Fig. 1. Entanglement between the two atoms (left plots) and between the two reservoirs (right 
plots) in the good cavity limit for two equal cavities with A = 0.1C and for two different dctunings: 
8 = Q/2 (solid, blue line) and 6 = 2Q (dashed, purple line). The initial state is such that s = 0.2. 



by the above expressions for the concurrences, we report in Fig. |T]) the dynamics of 
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atom-atom entanglement, Ci.2(i), and reservoir-reservoir entanglement, C riiT - 2 (t), 
for two different detunings and for equal cavities, whose parameters are very close 
to the ideal limit (the photon escape rate is chosen to be ten times smaller the 
resonant Rabi frequency). Various features should be noted in these plots; let us 
concentrate on the case with smaller detuning first (8 = 0/2, solid lines). Under 
such conditions, C\^ shows dark periods; that is, time intervals following an ESD, 
during which the entanglement stays equal to zero before abruptly reviving with 
a periodicity essentially dictated by the ideal Rabi time, 2ir/R. During these dark 
intervals, C ri:r - 2 attains its maxima, and the two concurrences oscillate exactly out 
of phase. On a time scale essentially set-up by A, both of the bipartite entangle- 
ment measures are seen to become close to their asymptotic values of C\ t 2 —> 0, 
and C riyr2 — > \/l — s 2 (the saturation value for the inequality reported above). 

One particular aspect that we want to consider here is the effect of the detun- 
ing. It was noticed in Ref. ^that, in the ideal Jaynes-Cummings limit, ESD is quite 
sensitive to the detunings. The same statement is clearly seen to hold for the reser- 
voir entanglement and for the ESB. Indeed, Fig. (pj shows that by increasing <5, 
apart from an obvious change in the oscillation period, the time evolution is qual- 
itatively modified as i) the ESD does not occur anymore (and, therefore, no dark 
intervals are found); ii) it takes a much longer time before the asymptotic values are 
approached, both for the atom-atom and for the reservoir-reservoir entanglement. 

This sensitivity to detunings, however, is seen to be gradually lost as the cavities 
become more and more leaky. As an example, Fig. @ reports the same plot as Fig. 
(U}, but for the case = A. It can be seen that, although the time required to 
almost reach the stationary values is again longer for larger detunings, the huge 
effect of 5 on the sudden death is not found anymore: the entanglement collapses 
to zero in both cases and only the ESD-time is affected by S. By further increasing 




Fig. 2. Same as Fig. |(T}, but with A = Cl. Here too, the solid, blue plots correspond to <5 = f2/2; 
while the purple, dashed ones correspond to the S = 2f2-case. 



the ratio A/fi (which implies going towards the bad cavity limit), the effect of 6 
becomes less and less noticeable. 

A similar behavior is found to occur when the two atoms arc differently de- 
tuned from their cavities. The effect is particularly enhanced when the detunings 
are opposite in sign (this is reminiscent of what has been found in Ref. ™). Indeed, 
by controlling the detuning S — S± = —62, one can qualitatively change the entan- 
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glement dynamics and explore regions in which ESD occurs, others in which dark 
periods or entanglement oscillations arc found, eventually reaching parameter re- 
gions in which a monotonic long-time decay takes place. As an example, in Fig. (|3j). 
we report four different behaviors of the atom-atom entanglement, obtained, for an 
intermediate coupling regime with the cavities, by just changing the detunings. 



l.Q- C 




Fig. 3. Atom-atom entanglement, Ci,2(£), evaluated in the case of equal cavities with Q = 2.5A 
and with an atomic initial condition set by s = 0.3. The four plots (starting from below) refer 
to symmetrically detuned atoms with 5\/Q. = — &2/Q = 0, 1,2,3 (solid, dashed, dotted and dot- 
dashed lines, respectively). 



4. Entanglement Generation 

In this section we describe some schemes for generating entanglement among one 
of the atoms and the two cavities, between the two cavity fields, and between one 
atom and a distant cavity. These schemes are based on the assumption that one can 
control the couplings between each atom and its surrounding cavity field as well 
as the atom-cavity interaction time. Both these requirements are experimentally 
feasible, e.g., in ion cavity QED systems, since the effective atom-cavity coupling is 
achieved via laser driven two-photon processes. The effective atom-cavity coupling 
in this case is given by g C R = gQ/A, with g the cavity coupling constant, f2 the Rabi 
frequency of the laser driving the Ra man transition, and A the laser and cavity de- 
tuning from an excited third level L^. In this case, g e e can be controlled by changing 
the laser intensity and/or detuning while the interaction time is determined by the 
length of the laser pulse. 
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4.1. Generation of W-states 

Inserting Eq. (j4]) into Eq. ([3]), we can explicitly write the total state in terms of its 
five components, as follows 



ItM) 



l.9i52>QS)|0 fcl fc2 

/Ci/C2 



-^[Si{t)£ 2 {t)\e 1 e 2 ) ®\Q kl Q k2 

£i(t)J 7 2 (t) \e1g2) ® |0 fel ,photon 2 ) 
Fi(t)£ 2 (t) \gie 2 ) (g> | photon x , 0fc 2 ) 
Fi (t)T 2 (t) \gig 2 ) <8> | photon^ photon 2 



(14) 



We note that, at a time t such that F\(jE) — 1 and £i(i) ~ 0, the state of the system 
approximates 



i^'(*)> — {y 



1-s 



I.92; 



|0 fcl fc2 > 



(15) 



A:i A'2 



+ e 



IC/3 



^2(*)|e 2 )(E> Iphoton^O^ 



(t)^2(t) |ff2> <8> I photon^ photon 2 ) | |. 9 i) , 



From this equation one sees immediately that for s = 1/3 and £ 2 (t) = -4= one 
obtains the W-state 



V2 



with 



\ilfw) = -^ (|010) + |100) + |001)) , 



|010) = \g 2 ) (g) |0 fel fe2 ) , 

&1&2 

|100) = |e 2 )(g)|photon 1 ,0 fe2 ), 

|001) = \g 2 ) |photon 1 ,photon 2 ) 



(16) 

(17) 

(18) 
(19) 



Assuming an intermediate coupling regime for the first cavity, e.g., Qi = Ai, and 
taking for simplicity Si = 0, we see from Eqs. (JS|)- ^ that J-i(t) ~ 0.99 at Ait = 3. 
The condition £ 2 (t) — -4= is then obtained for 2 — 0.4A 2 , with Ai = A 2 and 5 2 = 0. 



4.2. Generation of Bell states 

W-states are characterized by the property that bipartite entanglement can be 
obtained after the measurement of any of the three parties. In our case, measuring 
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the atomic state of the second atom in \g2) generates a maximally entangled state 
of the cavity fields of the type 



fel fe2 ) + Iphoton^photoi^) . (20) 



The generation of a W-state followed by a measurement of the second atom in its 
ground state is similar to entanglement swapping, since it amounts at transferring 
entanglement from the atoms pairs to two fields inside two distinct non-interacting 
cavities. 

A Bell state of the atom-cavity system of the second atom can be generated 
by measuring the presence of a photon in the other cavity. In this case the state 
generated is 

IVW 2 ) = A= ( |e 2 fc2 ) + \g 2 , P hoton 2 ) J . (21) 

Most intriguingly, finally, a measurement of the cavity 2 in the vacuum creates 
a maximally entangles state of the atom 2 with the cavity 1, 

\^a 2ri ) = -^ f (g) | ff2 fcl ) + |e 2> photon!) J . (22) 

Finally we note that a similar W-state can be obtained by requiring that J-" 2 (t) ~ 
1 and £2(0 — 0. The results are exactly the same as the ones discussed above 
provided that one exchanges atom 1 with atom 2. 

5. Summary and conclusions 

We have discussed various features related to the transfer of quantum correlation 
due to the interaction of a pair of atoms (qubit) with two independent bosonic 
reservoirs, having non-markovian features. In particular, we have taken a Lorcntzian 
spectral density for each reservoir, in ord er to mimic the nearly resonant coupling of 
initially entangled atoms to leaky cavities^^. In such a case, the Hamiltonian model 
is amenable to an exact solution, so that the time evolution of various bipartite 
entanglement measures can be easily evaluated. In particular, we focussed on the 
atom-atom and reservoir-reservoir entanglement to show that: 1) the initial atomic 
entanglement is found to be transferred to the reservoirs; 2) the dynamics of the 
entanglement is strongly sensitive to the detuning of the atoms from their cavities, 
so that the detuning can be used as a knob to qualitatively change the character 
of the time evolution of the quantum correlations. These conclusions hold true in 
any atom-cavity coupling regime, ranging from ideal cavity (strong coupling regime) 
with a monochromatic spectral density, to a markovian environment with an almost 
structure- less environmental spectral density (bad cavity or weak coupling limit). 
Although the sensitivity to the detuning diminishes by going towards the markovian 
limit, we have found that an enhanced sensitivity can be obtained by independently 
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controlling the two atomic dctunings; and, in particular, that the entanglement 
time-evolution can be substantially altered if the two atoms are symmetrically out 
of resonance with respect to their cavities. 

We have also discussed a protocol to dctcrministically generate a W-state with 
one atom and two cavity fields, which can subsequently be used as a starting point 
for probabilistic schemes aimed at generating maximally entangled states of one 
atom and one, distant cavity field, and of the two cavity fields. It is worth mentioning 
that these schemes arc designed by taking into account the cavity losses and that 
they work despite their presence. 
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